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Abstract 

In this work, we evaluate the energy spectra of baryons which consist of two heavy and 
one light quarks in the MIT bag model. The two heavy quarks constitute a heavy scalar or 
axial vector diquark. Concretely, we calculate the spectra of \q(QQ') >i/2 an d \q(QQ') >3/2 
where Q and Q' stand for b and/or c quarks. Especially, for \q(bc) >i/2 there can be a mixing 
between \q(bc)a >i/2 and \q{bc)\ >i/2 where the subscripts and 1 refer to the spin state 
of the diquark (be), the mixing is not calculable in the framework of quantum mechanics 
(QM) as the potential model is employed, but can be evaluated by the quantum field theory 
(QFT). Our numerical results indicate that the mixing is sieable 



1 Introduction 

At present, the non-perturbative QCD which dominates the low energy physics phenomena, is 
not fully understood yet, a systematic and reliable way for evaluating the non-perturbative QCD 
effects on such as the hadron spectra and hadronic matrix elements is lacking. Fortunately, 
however, for the heavy flavor mesons or baryons which at least contain one b or c quarks 
(antiquarks) the situation becomes simpler due to an extra SUf(2) <S> SU S (2) symmetry pp. The 
studies in this field provide us with valuable information about the QCD interaction and its low 
energy behavior. 

Among all the interesting subjects, the hadron spectra would be the first focus of attention. 
The spectra of the J/ ip and T families have been thoroughly investigated in different theoretical 
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approaches. Commonly, the spectra are evaluated in the potential model inspired by QCD, where 
the QCD Coulomb-type potential is directly derived from the one-gluon-exchange mechanism, 
and the confinement term originating from the non-perturbative QCD must be introduced by 
hand[2]- For the heavy quarkonia where only heavy quark flavors are involved, the potential 
model definitely sets a good theoretical framework for describing such systems where relativistic 
effects are small compared to the mass scale. An alternative model, the bag model can also 
provide a reasonable confinement for quarks. In fact, for light hadrons, especially light baryons, 
the bag model may be a better framework for describing their static behaviors. The MIT bag 
model has some advantages (Sj. First, even though it does not hold a translational invariance, 
the quarks inside the hadron bag obey the relativistic Dirac equation and moreover, they can 
be described in the Quantum Field Theory (QFT) framework, namely there exist creation and 
annihilation operators for the constituents of the hadron. The latter property is very important 
for this work, that we can calculate a mixing between \q(bc)o >\/2 an d \q{bc)\ >i/2 (the notations 
will be explained below) states and it is impossible in the Quantum Mechanics (QM) framework. 

Another interesting subject is if the diquark structure which consists of two quarks and 
resides in a color-anti-triplet 3, exists in baryons. Its existence, in fact, is till in dispute. For 
light diquark which is composed of two light quarks, the relativistic effects are serious and 
the bound state should be loose. By contraries, two heavy quarks (b and c) can constitute a 
stable bound state of 3, namely a diquark which serves as a source of static color field 4 . As a 
matter of fact, the un-penetrable bag boundary which provides the confinement conditions to the 
constituents of the hadron, is due to the long-distance non-perturbative QCD effects, to evaluate 
the spectra, one needs to include the short-distance interaction between the constituents and 
it can be calculated in the framework of the perturbative QCD. In this work, we are going to 
evaluate the spectra of baryons which contains two heavy quarks (b and/or c) and a light quark 
and take the light-quark-heavy-diqaurk picture which obviously is reasonable for the case of 
concern. 

For evaluating the hadron spectra, the traditional method is the potential model. For 
baryons, the quark-diquark picture can reduce the three-body problem into a two-body problem 
and leads to a normal Schrodinger equation. Solving the Schrodinger equation, one can get the 
binding energy of quark and diquark[Sl EJ. In recent years, remarkable progresses have been 
made along this direction. The authors of refs. jjj have carefully studied the short-distance and 
long-distance effects, then derived a modified potential and obtained the spectroscopy of the 
baryons which contain two heavy quarks by using the non-relastivistic Schrodinger equation. 
Meantime, in Ebert et al.'s new work, the light quark is treated as a fully relativistic object and 
the potential is somehow different from that in their earlier work ^8 . In their works, not only 
the ground states of such baryons are obtained, but also the excited states are evaluated. 

However, the potential model has two obvious drawbacks. First, even though the diquark is 
heavy, the constituent quark mass of the light quark is still comparable to the linear momentum 
which is of order of Aqcd- Thus the reduced mass is not large and the relativistic effects are still 
significant. Secondly, working in the framework of QM, it is impossible to estimate the mixing 
of \q(bc)o >i/2 and \q(bc)i >xii where the subscripts and 1 of the (bc)-diquark denote the total 
spin of the subsystem, i.e. the (bc)-diquark (we only consider the ground state of I = 0). The 
reason is that there are no creation and annihilation operators in the traditional QM framework, 
so the transition (bc)i + q — > (bc)o + q, i.e. A + q — » S + q where the notation A and S refer 
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to the axial-vector and scalar diquarks respectively, is forbidden, even though the transition is 
calculable in QFT. 

On other side, the MIT bag model does not suffer from the two drawbacks. In this picture, 
since the diquark is heavy, it hardly moves so that can be supposed to sit at the center of 
the bag, whereas the light quark is freely moving in the bag and its equation of motion is the 
relativistic Dirac equation with a certain boundary conditionjH] and both the quark and diquark 
are quantized in the QFT. Thus the relativistic effects are automatically included. Secondly, 
one can deal with a possible conversion of the constituents in the bag in terms of QFT, namely 
one can let a constituent be created or annihilated, thus the transition A + q ^ S + q is allowed 
and the corresponding mixing of \q(bc)o >i/ 2 and \q(bc)\ >\/2 is calculable. 

Usually the bag model is not very applicable to the light mesons because the spherical 
boundary is not a good approximation for the two-body system. Even though the quark-diquark 
structure is a two-body system, the aforementioned problem does not exist because the diquark 
is much heavier than the light quark. The picture is in analog to the solar system or an atom 
where only one valence electron around the heavy nucleus, and spherical boundary would be a 
reasonable choice. 

In this work, following the literature [3], we treat the short-distance QCD interaction be- 
tween the light quark and heavy diquark perturbatively. Since the interaction energy Ei n t(R) 
is not diagonal for \q{bc)\ >\/2 an d \q(bc)o >\/2, we may diagonalize the matrix to obtain the 
eigenvalues and eigenfunctions which would be the masses of the baryons with flavor q(bc) and 
spin 1/2. Moreover, for the other baryons \q(bb)\ >i/2(3/2) \l( cc )i >i/2(3/2) \l(bc)i >3/2> the 
diquark must be an axial vector due to the Pauli principle^- 

The paper is organized as follows, after the introduction, we derive all the formulation of 
Ei n t(R) and Mb in Sec. II, then in Sec. Ill we present the numerical results and all concerned 
parameters, finally the last section is devoted to discussions. 



2 Formulation 

1. A brief review of the MIT bag model 

The wavefunction of a light quark in the MIT bag obeys the Dirac equation for free fermion 
and a boundary condition which forbids the quark current to penetrate the bag boundary. It 
has a formjnj 

9(r ' tj " [ -(^) 1/2 3i(xr/R)cr-rU, J W 



with 

N-\x) = R^U) M "- 1/R)+ , m/R , (2) 

uj(u — m) 

where ji is a spherical Bessel function, U is a two-component Pauli spinor, the eigen-energy is 

[X 2 + (mit!) 2 ] 1 / 2 

uj{m,R) = , (3) 
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and the eigenvalue x satisfies an equation 

X 

taIl(X) = l- mj R-[ X 2 + (mjR) 2]l/2- W 

In our picture, the two heavy quarks constitute a diquark which is a boson-like bound state of 
color 3 and because it is heavy, it hardly moves, but sits at the center of the bag. Its wavefunction 
can be written asjlUj 

N Ar£ 

ijj(r) = — e 2 for the scalar diquark; (5) 
Air 

N at£ 

Vv( r ) = ~^~ e 2 V/i for the axial vector diquark, (6) 

where 7]^ is the polarization vector which is normalized as rj 2 = — 1, A is a parameter and will 
be discussed later in the text. 

2. Formulation for the baryon spectra. 

In the CM frame of the baryon, the total mass of the baryon can be written as 

M B = M D +u; + E int (R) + ^TTR 3 B-^ (7) 

where Ei n t{R) is the interaction energy between the diquark and light quark, B is the bag 
constant and z is a constant for the zero-point energy, u and Mp are the eigen-energy of the 
free light quark and the mass of the heavy diquark which are given in the literature [Sj- In this 
work, following the standard procedure 13]. we are going to calculate the interaction energy 
Eint{R)- Generally, the interaction hamiltonian in the bag can be expressed as 

H D , D = -|| 5s 2 / g(x) 7 ^(x)^ t ^-(y) < ^)£B > *(y)d 3 xd 3 y. (8) 

It is noted that without mixing, i.e. Hd'd is diagonal, E,^ = H^id, however, if there is non- 
diagonal H d'd, Ei n t is the eigenvalues of the hamiltonian matrix. The expectation value of the 
Casimir operator < 0|Af AflO >= —16/3, the strong coupling g 2 s = 47ra s , g(x) is the wavefunction 
of the free light quark, < D'\J U \D > is the effective vertex for DD'g and D^ u is the gluon 
propagator in the Coulomb gaugefT!?]. The form of such a propagator in the configuration space 
reads 

G(r ' ro) = 4^F^ + S i + (i + i)k ^ (flMfl)1 - (9) 

In this work, ro is small, thus the main contribution comes from the I = component. Then, 
the expression can be simplified as: 
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It is noted that the term is due to the mirror charge effect [T2*]. in fact, it indeed corresponds to 
the zero-point energy in the bag as Jaffe et al. suggested [H]. Actually, k is related to the vacuum 
property and still serves as a free parameter that cannot be determined from any underlying 
theory yet. We choose k = 0.87 to fit the most recent lattice result about the spectrum of 
(ccq) l+ . 

2 

3. The D'Dg effective vertex 

If we only consider the ground states of the diquark, namely the two heavy quarks are in 
/ = color-anti-triplet state, the diquark can be either a scalar (denoted as S) with s = or an 
axial vector (denoted as A) of s = 1. The effective vertices can be derived by the quantum field 
theory under the heavy quark limit |13|. The effective vertex for SS'g is 

< S'\J U \S >= ^MM'(f lV ' u + f 2 v u ) (11) 

and the AA'g effective vertex is of the form 

<A'\J U \A> = VMM'[f 3 ( v ■ ri'*)v' u + / 4 (r/* • rf)v v + h{r) ■ v')(rf* ■ v)v' u 

+ h(v ■ v'Kv'* ■ v)v v + fW^iv ■ v 1 ) + h{rf* ■ v)rf] (12) 

where v, v' , rj^ and rj^ are the four-velocities and polarization vectors of D and D' respectively, 
/j's are the form factors at the vertices. As we did in our previous work jH] where the effective 
potential model was employed, for convenience of calculation, we would write the polarization 
into a spin-operator which acts on the wavefunction of the axial-vector diquark as 

fyi = -/=(£■ s>s), (13) 

where higher order relativistic corrections proportional to -^rr are neglected and the factor l/y/2 

is a normalization factor because < s 2 >= s(s + 1) = 2. It is worth noticing that Sj and Sj do 
not commute with each other, thus one must be careful about their order in deriving formula. 
The effective vertex for ASg is written as 

<A'\r\S> = ^[MM'\h\fi* V + M?/* • v)v' v 

+ Mr/* ■ v)v" + fuie^v'fv'^]. (14) 

Here we must stress that for convenience of calculation, we turn rj^ into the quantum spin- 
operator as evaluating the interacting energy for \q(QQ') >, but as pointed above, we cannot 
calculate the mixing between \q{bc)\ >\/2 and \q(bc)o >\i% in QM, but need to carry out the 
derivation in QFT instead. Then, we have to keep the polarization r/„ in a 4-vector form as 

^ = =F^(0,l,±i,0) and 77° = (0,0,0,1) 

. Because the diquark is very heavy and hardly moves, \(3\ <C 1, one can use the polarization 
vector in the reference frame where the diquark is at rest. 
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In the heavy quark limit, we have 

h = /2 = h — h = ~h = —fi = /i4 = l 

hi = hi = hs = 0, (15) 

and for very heavy diquark the above approximation holds 13 . 
4. The interaction energy 

To calculate the interaction energy, the basic formula is eq.(J8J). For the heavy diquark, 

2 

|p <C Mjj, thus the relativistic corrections proportional to and higher than -fe- can be safely 
ignored. 

Based on the approximations, one can easily obtain the interaction energies. 

For the baryon where the diquark is a scalar of 3, the interaction energy between the 
light quark and the scalar diquark is corresponding to the transition matrix element of < 
q, S\H e jf\q, S >, and can be expressed as 

Hss = -yy / /g(x)7^(x)^^*(y)[/ 1 ^ + / 2 ^]^(y)d 3 xd 3 y 

= & /\ 2 (^)+^)]-d 3 z, (16) 
JO K K r x 

where x is the spatial coordinate of the light quark, y is that of the heavy diquark, r = x — y 
is the relative coordinate of the quark and diquark. 

For the baryon where the diquark is an axial vector, the matrix element < q, A\H e ff\q, A > 
is written as formula of H AA i : 

h aa\ = -yf / /?( x )7' 1 9W^**(y)[/3(')-V>"'+/4(f-'))/ 

+ MV ■ v')(v'* ■ v)v' v + U V ■ v')(r,'* ■ v)v" + f rV '*(v ■ ^0 

Jo R R J |r x — r y | R 

+ 9 2 s ^f f R U^h^(^)d 3 x /(SxqH-^— + A)|*(y)| 2 A(17) 
M Jo R R J |r x — r y | R 

where C a p refers to the spin projections of the quarks in the baryon. The transitional momentum 
q will change into the form — iV in configurational space acting on the relative position, and we 
then get: 

where q is the exchanged momentum between the quark and diquark, and in the configuration 
space of the bag, it is an operator acting only on the relative coordinate r as — zV r and 

• q d 3 q ir 



|q| 2 (2vr) 3 r 3 ' 



(19) 



6 



We finally obtain an integral which must be carried out numerically if we take the form of ^(y)) 2 
which is treated as 5 function into our consider, we will finally read: 

2 N' f R . , x r x \ . / X r x \ i - . r (cr - S) 

+ 9 S — L M—)Ji(—) < sim lS2 m 2 |[— — 



M Jo R R 
- {a ' rx) f' rx) }\sim lS2 m 2 > d 3 x (20) 

r x 

For the baryon of spin 3/2, which is composed of the light quark of spin 1/2 and an axial 
vector diquark of spin 1, the interaction energy is 

H AAi = 9 ; w /;t^)+^)]AA 

- '""f" 1 !^ .M>A (21) 

For a mixing hamiltonian which originates from the transition S+q — > -A +(7, the non-diagonal 
interaction element is 

= | |g(x) 7 ^(x) J D M ^*(y)(ie^r ? ^> (7 )^(y)d 3 x a ! 3 y (22) 
After a straightforward integration, we obtain 

- (g ' Fx) f • rx) ]| ai mi>d 3 x (23) 

It is noted that here 77^ remains as a four- vector, and it reduces into a three- vector rj which is 
not an operator. 

With the mixing between \q(bc)o >i/2 and \q(bc)i >i/2, one can write the real eigenstates of 
(qbc) 1/2 and (qbc)^^ as 

l|.5>=c|i.i,o,o> + fl(- 1 /I|i,i,i,o> + y||i,-i,i,i», (24) 

where C and D are the coefficients to be determined. The interaction energy is 

< > = < I,I,o,o|F|i,i,o,o > +J D 2 (i < i,I,i, |tf|i i 1,0 > 



211 , ,i 1 2V2 11 , 1 1 

+ - < -, — -, — .1,1 > — — < -,-,i,oLff -, — ,1,1 > 

3 2 2 1 '2 2 3 2 2 1 '2 2 
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/2 1 1 , r ,1 1 
+ V3 < -,--,l,l|i?|-,-,0,0>), (25) 

which contains both diagonal and non-diagonal interaction elements. It is easy to set it into a 
matrix form as 

h= [h A8 h aa )> (26) 

and the corresponding Schrodinger equation is 

Diagonalizing the matrix, one can solve C and D which would determine the fraction of |g(bc)o >i/2 
and \q(bc)i >\/2 in the eigenstates. 



3 The numerical results 

In this work, we have the input parameters as 
m u = 171,1 = m q pa 0, m s = 0.279 GeV, 
M cc = 3.26 GeV, M bb = 9.79 GeV, M bc = 6.52 GeV 0. 
q s =0.23, B = (0.145GeV) 4 . 

As argued above, the zero-point energy is not considered. Minimizing the expression of Mg 
in eq.(l) with respect to the bag radius R, we obtain an equation and then determine the R- value. 
Substituting the obtained R-value into the expressions, we achieve the following table. In the 
table, we list our numerical results for the baryons which contain two heavy quarks, meanwhile 
the results obtained in other approaches are presented in the table for a clear comparison. 
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N otation 


content 


J p 


M B 


M B 


M B 


M B 


Mr 


Mr 


Mr 








(our results) 


[5] 


[7] 


[8] 


[14] 


[15] 


[16] 


1 — 'CC 


(cc)a 


1 + 

2 


3.55 


3.66 


3.48 


3.620 


3.55 


3.66 


3.61 


1— 'cc 


(cc)q 


3 + 
2 


3.59 


3.81 


3.61 


3.727 


3.64 


3.74 


3.68 




(cc)s 


1 + 


3.73 


3.76 


3.59 


3.778 


3.66 


3.74 


3.71 


IT 

J "cc 


(cc)s 


2 


3.77 


3.89 


3.73 


3.872 


3.73 


3.82 


3.76 


1 — '00 




1 + 

2 


10.10 


10.23 


10.09 


10.202 




10.34 




"* 
^bb 


(bb)q 

V / 1 


3 + 
2 


10.11 


10.28 


10.11 


10.237 




10.37 




Qbb 


(bb)s 


1 + 

2 


10.28 


10.32 


10.21 


10.359 




10.37 




O* 
ll bb 


(bb)s 


3 + 
2 


10.29 


10.36 


10.26 


10.389 




10.40 




^cb 


(cb)q 


1 + 

2 


6.80 


6.95 


6.82 


6.933 




7.04 




^cb 


(cb)q 


1 + 

2 


6.87 


7.00 


6.85 


6.963 




6.99 




"c6 


(cb)q 


3 + 
2 


6.85 


7.02 


6.90 


6.980 




7.06 






(cb)s 


1 + 

2 


6.98 


7.05 


6.93 


7.088 




7.09 




Kb 


(cb)s 


1 + 

2 


7.05 


7.09 


6.97 


7.116 




7.06 




Kb 


(cb)s 


3 + 
2 


7.02 


7.11 


7.00 


7.130 




7.12 





Table 1. The baryon spectra 



4 Conclusion and Discussion 

In this work, we evaluate the spectra of baryons which contains two heavy quarks in the MIT bag 
model. It is an approach parallel to the potential model which is widely adopted for studying 
heavy hadrons. One can notice from the table given in the text that the numerical results in the 
two approaches are consistent with each other by the order of magnitude, but there are obvious 
distinction in numbers. 

The confinement in both the potential model and the bag model is artificially introduced 
which may reflect the non-perturbative QCD behavior in certain ways, but since none of them 
are derived from the first principle, one cannot expect them to be precise and it is reasonable that 
they result in different numbers which are related to the physics pictures and phenomenological 
parameters. As a matter of fact, in general the parameters are obtained by fitting data and 
while calculating the spectra of the lowest states of heavy quarkonia, various potential forms 
and sets of parameters can meet data. The MIT bag model is an alternative model and this 
work may be complimentary to the potential model. 

No matter in the potential model or the bag model, there is a zero-point energy problem, 
which manifest the vacuum property of QCD and is not calculable at present. The zero-point 
energy would determine the exact positions of each baryon and its excited states, but dies not 
influence the relative distances between the states. In the bag model, situation is a bit different, 
because Eq oc 1/R, when to obtain Mb, we differentiate Mr with respect to R and the minimum 
determines the R- value. Since the energy of free light quark is \/x 2 + (rnr)' 2 /R has the same 
form of Eq, thus its contribution can be attributed to the little shift of the quark mass. The gaps 
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between various states are not affected. In the future when data of the spectra are accumulated, 
one can come back to make more accurate adjustment of all the phenomenological parameters. 
In this work, we ignore the zero-point energy as argued above. 

The advantages of using the bag model to evaluate the spectra are two- folds. First the 
light quark obeys the relativistic Dirac equation, secondly, one can use the QFT to evaluate the 
possible mixing between \q(bc)o >\/2 and \q{bc)\ >\/2 which is impossible in the regular QM 
framework. In our case, when we calculate the relativistic corrections to the interaction between 
quark and diquark, we only keep the terms up to order p/M, which results in splitting of Bi/ 2 
and B 3 / 2 of the same flavor combination. Because the diquark is very heavy and M S> |p|, this 
approximation would not bring up significant changes. 

As a conclusion, we evaluate the spectra of heavy baryons containing two heavy quarks 
in the MIT bag model. The results are qualitatively consistent with that obtained in the 
potential model, but numerically differ by a few percents. Moreover, we obtain a relatively large 
mixing between \q(bc)o >\/2 and \q{bc)\ >\/2 which will be tested in the future measurements. 
Since there are no data available so far, we cannot fix a few parameters such as the zero-point 
energy, and one can definitely expect some deviations from the real values. Once the data are 
accumulated in the future, we can fin-tune the model and parameters. Definitely, the data will 
provide us with valuable information about the model and parameters. 

Even though the present B-factories BaBar and Belle cannot produce such baryons, TEVA- 
TRON and LHC which will begin running in 2007, may measure the spectra of such baryons, 
especially the long-expected Next-Linear-Collider (NLC) will offer us an idealistic place for 
studying such hadron spectra and our knowledge on the hadron structure can be greatly en- 
riched. 
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of China. 
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